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We report on recent calculations of the differential cross section for top-quark pair production at hadron 
colliders. The results are differential with respect to the top-pair invariant mass and to the partonic scattering 
angle. In these calculations, which were carried out by employing soft-collinear effective theory techniques, we 
resummed threshold logarithms up to next-to-next-to-leading logarithmic order. Starting from the differential 
cross section, it is possible to obtain theoretical predictions for the invariant-mass distribution and the total cross 
section. We summarize here our results for these observables, and we compare them with the results obtained 
from different calculational methods. 



1. INTRODUCTION 

The study of top-quark properties and in par- 
ticular the measurements related to the produc- 
tion of top-quark pairs play a crucial role in the 
physics programs at the Tevatron and at the 
LHC. The large number of top-quark pair events 
which will be detected at the LHC will turn the 
study of these observables into precision physics. 
Precise measurements require equally precise the- 
oretical predictions for the quantities which are 
measured. Currently, complete next-to-leading 
order (NLO) QCD calculations are available for 
the total top-quark pair-production cross section 
[1-4] , as well as for differential distributions [5-7] . 
The predictions obtained from these calculations 
show theoretical uncertainties larger than 10 % 
for both Tevatron and LHC center of mass en- 
ergies. The calculation of the NNLO corrections 
for the production of top-quark pairs would be 
extremely useful. However, at present the two- 
loop virtual corrections which are needed to ob- 
tain complete NNLO calculations are known only 
in part [8-10]. It is possible to improve the cur- 
rent NLO predictions by employing threshold re- 
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summation methods. For the total top-quark pair 
production cross section, complete next-to-next- 
to-leading logarithmic (NNLL) expression were 
obtained very recently [11-13]. For what con- 
cerns the distributions differential with respect 
to the top-quark pair invariant mass and scat- 
tering angle, a NNLL resummation scheme was 
presented in [14], while approximate NNLO for- 
mulas were first published in [15]. Both calcula- 
tions were made possible by the evaluation of the 
needed two-loop anomalous-dimensions matrices 
[16,17]. In this work we summarize and comment 
the results reported in [14,15]. 

2. DIFFERENTIAL CROSS SECTION 



We consider the process 
AW) A 2 (P 2 ) -> t(p 3 ) i(p 4 ) X( PX ) , 



(1) 



where N\ , A 2 indicate either two protons (LHC 
case) or a proton and an antiproton (Tevatron 
case), while X indicates an inclusive hadronic fi- 
nal state. The doubly-differential cross section 
for the process in Eq. (1) can be written as 
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where M is the pair invariant mass, 9 is the scat- 
tering angle in the partonic frame, m t is the top 
quark mass, ff+j is the luminosity for the partonic 
initial state ij, s (s) is the hadronic (partonic) 
center of mass energy squared, and finally 
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The functions dj in (2) are the hard-scattering 
kernels which describe the partonic process. We 
define the partonic threshold as the region of the 
phase space in which s —> M 2 (z — > 1). The par- 
tonic threshold does not coincide with the produc- 
tion threshold region, s — > 4m 2 , which is usually 
considered in relation to the calculation of the to- 
tal pair-production cross section. At the partonic 
threshold there is no phase space available for the 
emission of additional hard radiation in the final 
state, therefore higher order corrections in that 
region are dominated by virtual corrections and 
by the emission of soft gluons. Partonic processes 
which are not initiated by a quark-antiquark pair 
or by two gluons are also suppressed at the par- 
tonic threshold. By using soft-collinear effective- 
theory methods, it is indeed possible to prove that 
in the z — > 1 limit the hard-scattering kernels can 
be written as [14] 
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where ij £ {qq, gg}; the hard functions H and the 
soft functions S are matrices in color space. (All 
of the objects in Eq. (4) also depend on the ar- 
guments m t ,cos#, and on the factorization scale 
fif.) The hard functions are related to virtual cor- 
rections, while the soft functions originate from 
the real emission of soft gluons. At n-th order 
in perturbation theory, the soft functions involve 
plus distributions of the form 
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In [15] we combined the explicit calculation of the 
hard and soft functions up to NLO with the infor- 
mation that can be extracted from the two-loop 



anomalous dimension matrices which appear in 
the rcnormalization group equations (RGE) sat- 
isfied by H and S [16,17]; in this way we ob- 
tained approximate NNLO formulas for the differ- 
ential cross section in (2). Such formulas include 
the analytic expression of the coefficients multi- 
plying the plus distributions in the NNLO hard- 
scattering kernels, the scale dependent terms mul- 
tiplying the delta function of argument (1 — z), 
and a specific set of subleading terms associated 
with the plus distributions 2 . In [14] we studied 
the resummation of the distributions in Eq. (5) 
up to NNLL directly in momentum space. (The 
same procedure was already successfully applied 
to other processes [18-23].) This was done by 
solving the RGE satisfied by H and S. The 
explicit formulas for the NNLL resummcd hard- 
scattering kernels can be found in [14], here we 
just want to remind the reader that they involve 
two extra scales on top of the factorization scale: 
These are the hard scale \ih and the soft scale \x s 
which characterize the hard and soft functions, 
respectively. At the moment of evaluating the 
resummed formulas numerically, /Lt s , (ih, and (if 
are chosen independently from one another. Both 
the approximate NNLO formulas and the NNLL 
resummcd formulas improve the fixed order calcu- 
lations of the hard-scattering kernels at NLO be- 
cause they include higher order corrections origi- 
nating from the partonic threshold region. Does 
the partonic threshold region provide a numeri- 
cally dominant contribution to the cross section 
in Eq. (2)? Since for phenomenologically interest- 
ing applications r < 0.3, at first sight this does 
not seem to be the case. However, if the luminosi- 
ties ff become small sufficiently fast for z r, 
the integrals in Eq. (2) will still be dominated by 
the region where z ~ 1. The latter effect goes 
under the name of dynamical threshold enhance- 
ment. By comparing the exact NLO invariant 
mass distribution with the approximate NLO in- 
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variant mass distributions obtained by retaining 
only the terms which are singular in the thresh- 
old region (and the sublcading terms associated 
to them), it is possible to show that a dynamical 
threshold enhancement docs take place in the pro- 
duction of top-quark pairs at the Tevatron and at 
the LHC [14,15]. 

3. PHENOMENOLOGY 

To obtain the best possible predictions, wc 
matched the NNLL resummed formulas to the 
NLO calculations for the invariant mass distri- 
bution and the total cross section. Also the ap- 
proximate NNLO corrections which we obtained 
are always added to the exact NLO results. In the 
NNLL+NLO calculations, the numerical values of 
the scales are chosen in order to minimize the de- 
pendence of the NNLL corrections on fj, s , fih, and 
fif, in particular fif and fin are set equal to M, 
while [i s is typically chosen to be in the interval 
[M/4, M/10] (see[14] for details). For the same 
reasons, when dealing with approximate NNLO 
corrections, we set fif = M. In general, once the 
NNLL or the approximate NNLO corrections are 
taken into account, the residual scale uncertainty 
on the top-quark pair observables is smaller than 
or comparable to the luminosity uncertainty. 

3.1. Invariant Mass Distribution 

The invariant mass distribution can be ob- 
tained by integrating the doubly differential dis- 
tribution of Eq. (2) in the range — 1 < cos < 1. 
The exact NLO results needed for the match- 
ing with the NNLL resummed and approximate 
NNLO corrections can be obtained by using the 
partonic Monte Carlo MCFM [26]. The results 
of the calculation of the invariant mass distribu- 
tions at the Tevatron are shown in Fig. 1. (Sim- 
ilar plots for the LHC case are available in [14].) 
All the bands in Fig. 1 have been obtained by us- 
ing MSTW2008 PDFs [24]. In the top panel, the 
LO band was obtained by employing LO PDFs, 
the NLO band by using NLO PDFs, etc. In the 
bottom panel, the NLL band was obtained by em- 
ploying NLO PDFs and the NNLL band by em- 
ploying NNLO PDFs. The width of the bands 
reflects the scale uncertainty of the distributions. 
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Figure 1. Top panel: Fixed-order predictions for 
the invariant mass spectrum at LO (light), NLO 
(darker), and approximate NNLO (dark bands) 
for the Tevatron. Bottom panel: Correspond- 
ing predictions at NLL (light) and NLO+NNLL 
(darker bands) in resummed perturbation theory. 

As expected, the width of the bands is signifi- 
cantly smaller in the lower panel, which shows 
the resummed results. The partonic threshold re- 
gion corrections accounted for in our formulas be- 
come more and more dominant as the value of M 
increases. The CDF collaboration released data 
for the invariant mass distribution in [25]. We 
compared the data with our NLO+NNLL calcu- 
lation and we found a good agreement, especially 
at large values of M (see Fig. 2). It is possible to 
rewrite the invariant mass distribution in terms 
of a distribution with respect to the variable (3 t 
in Eq. (3). The /^-distributions show clear peaks 
for 0t ~ 0.6 — 0.8 at both the Tevatron and at 
the LHC. We observe that j3 t is always smaller 
than f3 = yl — 4m^/s. Therefore, values of the 
total cross section which are obtained by integrat- 
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Figure 2. Comparison of the RG-improved pre- 
dictions for the invariant mass spectrum with 
CDF data [25]. 



ing the P — > limit of the partonic cross section 
should be taken with a certain degree of caution, 
as we will illustrate further in the Section 3.2. 

3.2. Total Cross Section 

The total top-quark pair-production cross sec- 
tion can be calculated by integrating the invariant 
mass distribution in the range 2rrit < M < i/s. 
An analytic NLO result for the total cross sec- 
tion can be found in [4]. In order to match our 
NNLL resummed results with the NLO cross sec- 
tion in [4] , we chose fit — 400 GeV, which is 
representative of the values of M at which the 
invariant mass distribution peaks. Our results 
for the total cross section are summarized in Ta- 
ble 1. In the second line of Table 1 we report 
the results obtained by retaining only the lead- 
ing singular terms in the threshold limit z — > 1 
in the expression of the NLO hard scattering ker- 
nels. By comparing those results with the ones 
for the exact NLO cross section (first line in Ta- 
ble 1) it is immediate to see that the neglected 
terms, which are sublcading in (1 — z), contribute 
only a few percent of the total cross section at 
NLO. This does not represent a proof of the fact 
that the threshold expansion works well also at 
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Figure 3. The 0(a s ) corrections to da/dfi at the 
Tevatron. Here /i/ = nit. 



higher orders in perturbation theory, however, it 
is quite reassuring. Our best predictions for the 
total cross section are obtained by matching the 
NLO calculations with the resummed formulas at 
NNLL and can be found in the last line of Table 1. 
The resummation has a relatively small effect on 
the cross section (10 — 15% enhancement), but 
it reduces significantly the scale dependence of 
the predictions. The enhancement of the cross 
section due to resummation can be mimicked in 
the NLO calculation by choosing a smaller fac- 
torization scale. The values obtained by choosing 
fif = m t are shown in Table 2. While the re- 
summed predictions are quite similar to the ones 
listed in Table 1, the NLO predictions in Table 2 
are significantly higher than the ones in Table 1. 

3.3. The Small /3-Expansion 

Recently, approximate NNLO calculations of 
the total partonic cross section which are valid 
in the production threshold limit (3 — > became 
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Table 1 

Results for the total cross section in pb. The first error refers to the perturbative uncertainties associated 
with scale variations, the second to PDF uncertainties. The factorization scale is fixed at [if = 400 GeV. 



H f = 400 GeV 


Tevatron 


LHC (7 TcV) 


LHC (14 TcV) 


t^NLO, leading 
^NLO + NNLL 


e 7Q+0. 79+0.33 
°- ' J -0.80-0.22 
c /(Q+0.78+0.33 
d -^ J -0.78-0.20 
O qn + o. 19 + 0.31 
u - ou -0. 19-0.23 


i oq+21+7 
— 19—7 

i34±i?±; 

149™ 


n a i +105+26 
' ol -101-27 
7fii +64+25 
' 0i -75-26 
OO1+40+24 
°^ i -42-31 



Table 2 

Results for the total cross section in pb. The first error refers to the perturbative uncertainties associated 
with scale variations, the second to PDF uncertainties. The factorization scale is fixed at \if = 173.1 GeV. 



fj, f = 173.1 GeV 


Tevatron 


LHC (7 TeV) 


LHC (14 TeV) 


Cnlo 

^NLO, leading 
^NLO+NNLL 


c 79+0.36+0.37 
u - ' z -0.76-0.24 
R ,49+0.42+0.35 
u '^ z -0.76-0.23 
c /fo+0.17+0.32 
°-^°-0. 21-0. 25 


159+*° + ? 


oon+107+31 
ooy -106-32 
oqr + 18+29 
ood -60-30 

jjl q + 50 + 30 

oio -36-35 



available [27,28,13]. These approximate formulas 
include logarithmic terms due to the emission of 
soft gluons and terms enhanced by Coulomb sin- 
gularities. By multiplying those expressions by 
the luminosity functions and subsequently inte- 
grating over 4toj < s < s it is possible to obtain 
values for the total cross section. The leading sin- 
gular contributions in the (3 — ¥ limit do not co- 
incide with those arising in the z — > 1 limit. How- 
ever, by expanding our approximate NNLO for- 
mulas for the scattering kernels, which are valid 
in the partonic threshold region, in the s — > 4m| 
limit and subsequently integrating over the scat- 
tering angle, it is possible to recover the results 
of [13] up to potential gluon terms. It is inter- 
esting to compare the predictions for the total 
cross section obtained by integrating the formu- 
las which are valid in the (3 -+ limit with the 
ones obtained by integrating the formulas which 
arc valid in the z — > 1 limit. This analysis can 
be carried out at NLO where exact results are 
available. We consider then the product of the 
luminosity and of the 0(a s ) corrections to the 
partonic cross section as a function of /?. Those 
quantities are shown in Figs. 3 (Tevatron case) 
and 4 (LHC case). In the figures, the solid black 
lines represent the exact NLO result, the solid 
red lines are our approximate NLO results (NLO 
leading), the dashed lines correspond to the ap- 



proximate NLO results in which we drop all of the 
sublcading terms (NLO leading PIM), and finally 
the dotted lines indicate the approximate NLO 
formulas in the (3 — > limit (NLO /3-exp, see [14] 
for further details). The contribution of each of 
these corrections to the total cross section is given 
by the area below the corresponding curve (the 
maximal value of /3 allowed is ~ 0.98 at the Teva- 
tron and ~ 0.99 at the LHC). As expected, all the 
curves become indistinguishable as (3 — > 0. How- 
ever, the NLO leading curves reproduce the shape 
of the exact results also for (3 3> 0, while the the 
NLO /3-exp curves do not. In particular, the NLO 
/3-exp overestimates the contribution to the cross 
section in the quark- ant iquark channel at both 
the Tevatron and the LHC. The NLO leading 
curves underestimate the corrections in the gluon 
channel; however, the good agreement between 
the area under the exact NLO curve and the area 
under the /3-exp curve at the LHC looks rather 
accidental. Finally, by comparing the NLO lead- 
ing curves to the NLO leading PIM curves, one 
sees that the subleading terms which are kept in 
the former sensibly improve the agreement with 
the exact NLO results. 
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Figure 4. The 0(a s ) corrections to da/d(3 at the 
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